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Abstract 

We develop a kind of fractional calculus and theory of relaxation and diffusion equa- 

t 

tions associated with operators in the time variable, of the form (J$n.\u)(t) = J k(t — 

o 

t)u(t) dr — k(t)u(0) where A; is a nonnegative locally integrable function. Our results are 
based on the theory of complete Bernstein functions. The solution of the Cauchy problem 
for the relaxation equation H>(k)U = —An, A > 0, proved to be (under some conditions 
upon k) continuous on [(0, oo) and completely monotone, appears in the description by 
Meerschaert, Nane, and Vellaisamy of the process N(E(t)) as a renewal process. Here 
N(t) is the Poisson process of intensity A, E(t) is an inverse subordinator. 
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1 INTRODUCTION 



The basic ingredient of the theory of fractional evolution equations (TJ [13], the Caputo-Dzhrba- 
shyan fractional derivative of order a G (0, 1), has the form 



t 

d 



dt 

where 



(B {k) u)(t) = - I k(t- r)u(r) dr - k(t)u(0) (1.1) 



o 



k(t) = — t>0. (1.2) 

1 (1 — a) 

A recent theory of evolution equations with distributed order derivatives (used for modeling 
ultraslow relaxation and diffusion processes; see jU El [191 ED] and references therein) is based 
on the operator (1.1) with 



Ht) = j ——-dp(a), t>0, (1.3) 
o 

where p is a Borel measure on [0, 1]; see [D il OH M, DSl OH EH M, EH]. 

It is natural to look at a general operator (1.1) and to ask the following question. Under 
what conditions upon a nonnegative function k G L 1 1 oc (lR + ) does the operator possess a 
right inverse (a kind of a fractional integral) and produce, as a kind of a fractional derivative, 
equations of evolution type? The latter means, in particular, that 

(A) The Cauchy problem 

(B( fc) u)(*) = -Xu(t), t>0; u(0) = l, (1.4) 

where A > 0, has a unique solution u\, infinitely differentiable for t > and completely 
monotone, that is (-l) n uf \t) > for all t > 0, n = 0, 1,2,.. .. 

(B) The Cauchy problem 

(B( fe )tu)(t, x) = Aw(t, x), i>0,i6l"; w(0, x) = w (x), (1.5) 

where wo is a bounded globally Holder continuous function, that is |u>o(£) — Wq(t))\ < — r]| 7 , 
< 7 < 1, for any £,r] G M n , has a unique bounded solution (the notion of a solution should 
be defined appropriately). Moreover, the equation in (1.5) possesses a fundamental solution of 
the Cauchy problem, a kernel with the property of a probability density. 

Note that the well-posedness of the Cauchy problem for equations with the operator D(^) 
has been established under much weaker assumptions than those needed for (A) and (B); see 

In the above special cases (A) and (B) are satisfied; see [3 [T3]. When k(t) has the form 
(1.2), the function u\ can be expressed via the Mittag-Leffler function E a : u\(t) = E a (—\t a ); 
see [13]. In the case (1.3), the asymptotic behavior of u\(t) is studied in [HI [16]; for the 
operator-theoretic meaning of the distributed order derivative and integral see [15] . 

From the point of view of mathematical physics, it is natural to expect the emergence 
of mathematical theories, in which the relaxation function u\(t) appears instead of e~ A *. So 
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far, the first developments of this kind are in the theory of stochastic processes, namely the 
renewal processes with slowly decaying distribution functions of waiting times; see [HJ [23] and 
references therein. 

In this paper we find a class of general operators (1.1) possessing the above evolution 
generating properties. This class is described in terms of analytic properties of the Laplace 
transform 

oo 

K{p) = J e- pt k(t)dt. (1.6) 
o 

We develop, in particular, a theory of the Cauchy problem (1.5). In contrast to the classical 
theory of parabolic equations and its analogs known for the cases (1.2) and (1.3), the main 
technical tools are not the contour integration and explicit estimates, but the theory of com- 
plete Bernstein functions [27]. Our solution of the Cauchy problem (1.4) leads to an analytic 
description of general renewal processes constructed (for a slightly more general situation) in 
[23J in terms of the random time change in the classical Poisson process determined by an 
inverse subordinator. 

In Section 2, we give a survey of the results we need about complete Bernstein functions 
and Stieltjes functions. In Section 3, we introduce and study an analogue, for the general 
framework, of the fractional integration operator. Section 4 is devoted to the problem (1.4), 
while Section 5 deals with the problem (1.5). 



2 Complete Bernstein Functions and Stieltjes Functions 

In this section we collect information we need about the classes of functions mentioned in the 
title. For the detailed exposition see [27] . 

A real-valued function / on (0, oo) is called a Bernstein function, if / G C°°, f(X) > for 
all A > 0, and 

> o for all n > 1, A > 0. 
Equivalently, a function / : (0, oo) — > R is a Bernstein function, if and only if 

oo 

/(A) = a + bX + J (1 - e~ xt ) fi(dt) (2.1) 
o 

where a, b > 0, and \i is a Borel measure on [0, oo), called the Levy measure, such that 

oo 

J mm(l, t) fx(dt) < oo. (2.2) 



o 



The triplet (a, b, /x) is determined by / uniquely. In particular, 



a = f(0+), b= lim (2.3) 

A— >oo A 
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A Bernstein function / is said to be a complete Bernstein function, if its Levy measure [i 
has a completely monotone density m{t) with respect to the Lebesgue measure, so that (2.1) 
takes the form 



/(A) = a + bX + I (1 - e~ xt ) m{t) dt (2.4) 

where, by (2.2), 



J min(l, t)m{t) dt 



< oo. 



Here the complete monotonicity means that m e C°°(0,oo), (-l) n m^\t) > 0, t > 0, for all 



n = 0,1,2, 



Another important class of functions is that of Stieltjes functions, that is of functions (p 
admitting the integral representation 

oo 

v (\) = ^ + b + J j^-f(dt) (2.5) 
o 

where a, b > 0, a is a Borel measure on [0, oo), such that 

oo 

/ (1 +t)~ 1 a(dt) < oo. 



Using the identity (A + t) 1 = J e ts e Xs ds we find from (2.5) that 

o 



(2.6) 



^(A) = - + 6+ / e- Xs g(s)ds (2.7) 



where 



oo 



9(s) = / e~ ts a(dt) (2.8) 



is a completely monotone function whose Laplace transform exists for any A > 0. 

We will denote the class of complete Bernstein functions by CBJ 7 , and the class of Stieltjes 
functions by S. The following characterization is proved in [27]. 

Proposition 1. Suppose that f is a nonnegative function on (0, oo). Then the following con- 
ditions are equivalent. 

(i) / G CBT. 

(ii) The function A >-)■ A _1 /(^) ^ n «S- 
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(iii) / has an analytic continuation to the upper half-plane H = {z G C : Imz > 0}, such that 
lmf(z) > for all and there exists the real limit 

/(0+) = lim /(A). (2.9) 

(0,oo)3A^0 

(iv) / has an analytic continuation to the cut complex plane C \ (— oo, 0], such that Imz ■ 
lmf(z) > 0, and there exists the real limit (2.9). 

(v) / has an analytic continuation to HI given by the expression 

oo 

f(z) = a + bz + J -i— cr(dt) (2.10) 
o 

where a, b > 0, and a is a Borel measure on (0, oo) satisfying (2.6). 



Note that the constants a,b are the same in both the representations (2.4) and (2.10). 
The density m(t) appearing in the integral representation (2.4) of a function / G CET and the 
measure a corresponding to the Stieltjes function <^(A) = A _1 /(A) are connected by the relation 

oo 

m(t) = J e~ ts sa(ds). (2.11) 
o 

The importance of complete Bernstein functions is caused by the following "nonlinear" 
properties [27], quite unusual and having significant applications. 

Proposition 2. (i) A function f ^ is a complete Bernstein function, if and only ifl/f is 
a Stieltjes function. 

(ii) Let f, fx, f 2 e CET , <p, (fi,(p 2 e S. Then fop e S, cpof e S, f x of 2 e CET , <piO(p 2 e CET , 
(A + /)" 1 eS for any A > 0. 

It follows from Propositions 1 (ii) and 2 (i) that ^ / G CEJ 7 , if and only if the function 
/*(A) = A//(A) belongs to CEJ- '. Let us write its representation similar to (2.4), 



/*(A) = a* + b*X + / (1 - e~ xt ) m*(t) dt. 

Then 



o 



a 



0, if a > 0, 

, if o = 0; (2-12) 

b+ J tm(t) dt 



(i 



'0, if 6 > 0, 

b*={^ , if 6 = 0; (2-13) 

a+ J m(t) dt 
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3 Fractional Derivative and Integral 

Throughout the paper (except the uniqueness theorem in Section 5) we make the following 
assumptions regarding the Laplace transform (1.6) of the function k. 

(*) The Laplace transform (1.6) exists for all p > 0. The function /C belongs to the Stieltjes 
class S, and 

/C(p) — > oo, as p — > 0; /C(p) — > 0, as p — > oo; (3.1) 

pK,(p) — )■ 0, as p — )■ 0; pJC(p) — >■ oo, as p — >• oo; (3.2) 

By Proposition 1, the function p i— >■ pK,(p) is a complete Bernstein function. It follows from 
(2.3), (3.1), and (3.2) that the integral representations (like (2.5) and (2.7)) of the function /C 
have the form 



oo 



Kip) = J ^rfidt) 





and 

oo 

JC( p ) = J e - ps g{s)ds (3.3) 
o 

where 

oo 

g{s) = J e~ ts a(dt), 
o 

and the measure a satisfies (2.6). For the function p h-> pJC(p) we have 

oo 

pK,{p) = J (1 - e- pt ) m(t) dt, 
o 

and the limit relations from (3.1) and (3.2) show that 

oo oo 

J m(t) dt = J tm(t) dt = oo. (3.4) 



It follows from the uniqueness theorem for the Laplace transform that g(s) = k(s), so that 
the assumptions (*) imply the representation 

oo 

k(s) = / e- ts a(dt), < s < oo, (3.5) 



so that k is completely monotone. 

For each fixed s > 1, the function 1 1— > (1 + t)e~ ts is monotone decreasing on [0, oo), and its 
value at the origin is 1. It follows from (3.5), (2.6), and the dominated convergence theorem 
that k(s) — > 0, s — > oo. 
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Note that the conditions (*) are obviously satisfied for the kernel (1.2); under some condi- 
tions upon the weight measure p, they are satisfied also for the case (1.3). On the other hand, 
given a function JC satisfying (*), one can restore k by the formula (3.5). As a simple example, 
consider the complete Bernstein function p i-> log(l +p^), < /3 < 1 (see Example 15.4.59 in 
|27J), and set JC(p) = p^ 1 log(l + p&). Then fC(p) ~ p 13 ' 1 , as p — > 0, and tC(p) ~ fip" 1 logp, as 
p — > oo, so that the conditions (3.1) and (3.2) are satisfied. The above asymptotic properties 
are different from those corresponding to the cases (1.2) and (1.3). 

By Proposition 2, the function p i-> belongs to the Stieltjes class. Using (2.12), 

p)C{p) 

(2.13), and (3.3) we find its representation similar to (3.3), that is 

e~ ps x(s) ds 



p)C(p) 



where x(s) is a completely monotone function, x(s) — >■ 0, as s — > 0. Just as in (3.5), we get 
the representation 

oo 

x{t) = J e~ xt ri(d\), < t < oo, (3.6) 
o 

, ^7](d\) 

where I < oo. 

o 1 + A 

Let us consider the convolution 



(k*x)(t)= k(t-T)K(r)dr. 



By the construction of x, the Laplace transform of k * x equals -, so that 

p 

(k*x)(t) = l. (3.7) 

In other words, k and x form a pair of Sonine kernels. Such kernels were searched for 
from the 19th century; see [26] for a survey. The connection between the complete Bernstein 
functions and Sonine kernels seems never noticed, but can be useful due to the availability of 
tables of complete Bernstein functions |27j . 

Let us study, under the assumptions (*), the generalized fractional differentiation operator 
Dm of the form (1.1) and the generalized fractional integration operator 



(W)(*) = J x(t-s)f(s)ds. (3.8) 
o 

The operator is defined on continuous functions u, such that k * u is almost everywhere 

differentiable, for example, on absolutely continuous functions u. The operator is defined 
on Li oc (0,oo). 

The following result extends a well-known property of the Caputo-Dzhrbashyan fractional 
derivative (see Lemma 2.21 and Lemma 2.22 in [TBI). 
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Theorem 1. (i) If f is a locally bounded measurable function on (0, oo), then 3^)I(k)f = /■ 
(ii) If a function u is absolutely continuous on [0, oo), then (I^D^-u)^) = u(t) — u(0). 

Proof, (i) Let v = \ k) f . By (3.6) and (3.8), 

t oo 

v(t) = J f(s)ds J e- x ^ V (d\), 



which implies the inequality 

oo 

\v(t)\ < cj j (1 -e- xt )r](d\), 0<t<T<oo, (3.9) 
o 

(here and below the letter C denotes various positive constants). By (3.9), 

1 oo 

\v(t)\<CtJ V (d\) + C J (l- e -*)5(^), 
o 1 

so that v(t) ->■ 0, as t ->■ +0. 

t 

By (3.7), we have k*v=k*x*f= 1 * f, so that (k * v)(t) = J /(r) dr is absolutely 

o 

continuous, and B>^)V = f. 

(ii) We find, using (3.7) twice, that 



t s 

(I( fc) B (fc) u)(t) = J x{t-s) ^-J k(s - t)u{t) dr 



k{9) 



ds — u(0) 



t-0 



J k(t - 9 - t)u{t) dr d9-u(0) = j t J x{9) J k(t - 9 - t)u{t) dr - «(0) 

J 

t 

d d f 

— —(x*k * u){t) - «(0) = — / u(t) dr = u{t) - «(0). ■ 



4 Relaxation Equation and Renewal Processes 

Let us consider the Cauchy problem (1.4). 

Theorem 2. Under the assumption (*), £/ie problem (1.4) has a unique solution u\(t), contin- 
uous on [0, oo), infinitely differentiate and completely monotone on (0,oo). 



8 



Proof. Applying formally the Laplace transform, we find for the image u\ of a solution the 
expression 

that is pu\(p) = $ o {pfC{p)) where $(z) = -. By Proposition l(v), $ is a complete 

2 -|- A 

Bernstein function. It follows from the assumption (*) and Proposition 2(h) that p 1— >■ pu\(p) 
is a complete Bernstein function. Therefore u\ is a Stieltjes function, so that by (2.7), 

00 

^ a +6+ [e^g(s)ds, 



plC(p) + A p 



where g is completely monotone, possesses the representation (2.8) with the measure satisfying 
(2.6), and by (2.3), (3.1), and (3.2), 

a = hm — . = 0, b = lrm — = 0. 

p-^+o pH,{p) + A p^oo pl^{p) + A 

Thus, we have found that 

OO 

K{p) 



pK(p) + A 



e ps g(s)ds, 



and the identity (4.1) will hold if we set u\ = g. In addition, it follows from (3.2) that 

£(p) 1 

-, p — >■ 00. 



pK{p) + A p 

Taking into account the monotonicity of g, we may apply the Karamata- Feller Tauberian the- 
orem [S] and find that U\(t) — > 1, as t — Y +0. 

Now it follows from the uniqueness theorem for the Laplace transform that u\ is a solution, 
in the strong sense, of the problem (1.4). ■ 

Remark. It is instructive to see what happens if IC(p) satisfies weaker conditions, for ex- 
ample, if p h-> pK,(p) is a Bernstein function, but not a complete Bernstein function. It is seen 
from comparing (2.1) and (2.4) that the simplest example of this kind is obtained if we take 
JC (p) =p- 1 (l~e- p ). Then 

00 f 

Mp) = / ut)e- pt dt, ko(t) = r * ° " * " 

J 0, if t > 1. 

^ 

It is easy to check that 

- w , ( w(t) - w(t - 1), if t > 1; 
^ (fco) ; 1 ; lu(t)-u(0), if <t < 1, 
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and a general solution of corresponding problem (1.4) has the form 




u(t) 



c = const. This function is not continuous at the origin. It is continuous at t = 1, if c = 1, but 
even in this case it is not differentiable at t = 1. 

A probabilistic interpretation of Theorem 2 can be given on the basis of the results by 
Meerschaert, Nane, and Vellaisamy [23]. Let D(t) be a subordinator (see [3 J for the necessary 
notions and results from the theory of Levy processes), 

E [e" sD M] = e"**« 

where the Laplace exponent \I/ is a Bernstein function having the representation 



oo 

*(s) = bs + J (1 - e- ST ) <&{di 





with the drift coefficient b > and the Levy measure $, such that either b > 0, or $((0, oo)) = 
oo, or both. 

The process D is strictly increasing, thus it possesses an inverse function 

E(t) = inf{r > : D(r) > t}. 

Let N(t) be the Poisson process with intensity A. It is shown in [23j that N(E(t)) is a 
renewal process with independent identically distributed waiting times J n , such that 

P[J n > t] = E [e~ xm ] . (4.2) 

The Laplace transform of the expression in (4.2) is as follows: 

oo 

/ e~ st E [e- XE ^] dt = .. . (4.3) 

J 1 J s(A + *(s)) 

o 

It is interesting to know analytic properties of the function (4.2), and we can find them for 
the case where = s/C(s) is a complete Bernstein function (so that K is a Stieltjes function) 
satisfying (3.1) and (3.2). In this case the right-hand side of (4.3) is exactly the Laplace 
transform of the function it\(t) described in Theorem 2. Therefore, under these assumptions, 
the function (4.2) is continuous on [0, oo) and completely monotone. Earlier such properties 
were known [23] for the cases of stable subordinators and their distributed order extensions. 
In the latter case they were obtained from a rather complicated explicit integral representation 
[T4] for u\(t), so that here we have given a simpler proof based on the theory of complete 
Bernstein functions. 
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5 Bm-Heat Equation 



Let us consider the Cauchy problem (1.5). Applying formally the Laplace transform in t to 
both sides of (1.5) we obtain the following equation for the Laplace transform w(p,x) of a 
solution of (1.5): 

p)C(p)w(p, x) — JC(p)wo(x) = Aw(p, x), p>0, x G M n . (5.1) 

A bounded function w(t,x) will be called a LT-solution of the problem (1.5), if w is con- 
tinuous in t on [0, oo) uniformly with respect to x G M. n , w(0,x) = w (x), while its Laplace 
transform w(p,x) is twice continuously differentiate in x, for each p > 0, and satisfies the 
equation (5.1). 

Theorem 3. Suppose that the assumption (*) is satisfied. There exists such a nonnegative 
function Z(t,x), t > 0, x G IR n , a; / 0, locally integrable in t and infinitely differentiable in 
i^O, that 

Z(t,x)dx = l, t>0, (5.2) 



and for any bounded globally Holder continuous function u> , the function 

w(t,x) = J Z(t,x-£)w (£)d£ (5.3) 

is a LT-solution of the Cauchy problem (1.5). 

Proof. Consider the function 

g(s,p) = IC{p)e- splc{p) , s > 0, p > 0. (5.4) 

Since p \— > pfC(p) is a Bernstein function, the function p i-> e~ spK: ^ is completely monotone 
(see conditions for the complete monotonicity in Chapter 13 of [8]). By Bernstein's theorem, 
for each s > 0, there exists such a probability measure [i s {dr) that 

oo 

e -spK( P ) = I e -pr^ dr y ( 5 _ 5 ) 



The family of measures is weakly continuous in s. 
Set 

t 

G(s,t) = J k(t - r)fx s (dr). 
o 

By (3.3) and (5.5), the Laplace transform in t of the function G coicides with the function 



g(s,p) = j e- pt G(s,t)dt. 



o 
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On the other hand, it is seen from (5.4) that 



oo 



g{s,p) ds = -, 
V 







so that 



which implies the equality 



oo oo 

J e~ pt dt J G(s,t) ds=-, 



J G{s,t)ds = 1. (5.6) 
o 

We define Z by the subordination equality 

oo 

Z(t, x) = J (47rs)- n/ V^G(s, t) ds, x^O. (5.7) 



The equality (5.2) follows from (5.6) and properties of the fundamental solution of the classical 
heat equation. It follows from (5.7) and (5.6) that Z(t,x) is infinitely differentiate in x ^ 0. 
By the definition of G, the Laplace transform in t of the function Z, 

oo oo 

Z(p,x) = J (4rrs)- n/2 e- ]J £g(s,p) ds = /C(p) / {^s)- n/2 e-^ e~ spK{v) ds 







exists for p > 0; by Fubini's theorem, this implies the local integrability of Z(t, x) in t. Moreover, 
using the identity 2.3.16.1 from [25] we find that 



Z(p,x) = (27r)-t| x |i-f/C(p)(p/C(p))2(f- 1 )ir § _ 1 (|x|v / p/C(p)), x ^ 0, (5.8) 

where K v is the McDonald function. Note that, by virtue of Proposition 2, JC(p) ^ for any 
p > 0. 

Let us consider the function (5.3) starting from its behavior as t — > 0. Using (5.2) we write 



w(t, x) - w (x) = J Z(t, x — £) [w (C) ~ wo(x)] d£, 

so that 

oo oo 

\w(t,x)-w (x)\ <C J s- n/2 G(s } t)ds J e~^r\z[ 1 dz = C l J s l/2 G{s,t) ds, 

E™ 

and in order to prove that w(t, x) — > w (x) (uniformly in x), as t — > 0, it suffices to show that 
the function 



a(t) = J s^ /2 G(s,t) ds 



o 
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tends to 0, as t — > 0. Consider its Laplace transform 

oo oo 

o(p)= / s' y/2 g(s,p)ds = K(p) [ s l/2 e- spK{p) ds = hp'^pKip)}- 1 ' 2 







where h > 0. 

Since the function A H- A 7 / 2 is a complete Bernstein function, we find from Proposition 2 
that p i — y [p/C(p)] 7 ^ 2 is a complete Bernstein function, the function p v-t [p/C(p)] _7//2 is a Stieltjes 
function possessing a representation of the form (2.7). Thus, 



[p/C(p)P /2 = ^ + d + J e- ps g^(s) ds (5.9) 



o 



where g 1 is a completely monotone function of the form (2.8) with the measure a satisfying 
(2.6), c,d > 0. 

To find the constants c and d in (5.9), we can consider the complete Bernstein function 



oo 

p{plC(p)]~ l/2 = c + dp + p J e- ps gi {s) ds 



and use the formulas (2.3). It follows from (3.1) and (3.2) that 

c = lim p[plC{p)}-^ 2 = lim p 1 -* \K,(p)}-^ 2 = 0, 
p— >+o p— >+o 

d = lim [pJC{p)}-^ /2 = 0. 

p— y oo 

We have found that 



a(p) = hp 1 J e ps g 1 (s)ds = h J e pt dt J g 7 (r) dr, 
o oo 



so that 



z 

a(t) = h J (? 7 (r) dr — > 0, as t — > 0. 



o 

In order to perform the Laplace transform of the both sides of (5.3), we need estimates of the 
function (5.8). It is known [2] that the function K v and all its derivatives decay exponentially 
at infinity. If v > is not an integer, then 

K v {z) = —?—[I- v {z)-I v {z)\ 
2 sin z/7T 

where I u (z) = z u (p u (z), and (p u is an entire function. Therefore near the origin 

\K2\z)\<C\z\-"- 1 , 1 = 0,1,2,.... (5.10) 
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If v is a nonnegative integer, then 



K v {z) = /„(*) log - + z-"P v {z) + Q v {z) 

where P u is a polynomial, Q v is an analytic function on a neighbourhood of the origin, so that 
the inequality (5.10) remains valid for all natural numbers v, and 



< 



c\\og Z \, a 1 = 0, 



C\z\~ l , if z > 1. 

As a result, it follows from (5.8), (5.10), and (5.11) that for each fixed p, j — 1, . . . , n, 



(5.11) 



& 



dx 



T Z(p,x) 



< C\x\ 



-n+2—l—a\x\ 



, / = 0,1,2,..., 



a > 0, if n > 3. If n — 2, then (5.12) remains valid for / > 1, while 



Z(p,x) 



< C|log|x||e- aW , n = 2. 



If n — 1, there is no singularity at x = 0: 



<Ce- a|a:| , n = l. 



Now the Laplace transform in (5.3) is legitimate, so that 



w 



ip, x) = J Z(p, x - 0^o(0 P > 0, 



(5.12) 



(5.13) 



(5.14) 



The estimates (5.12)-(5.14) justify a single differentiation in spatial variables: 
dw(p,x) f d 



dxj 



dxj 



Z(p,x - £)w (£)d£, j = l,...,n. 



If n — 1, also the direct second differentiation is possible, thus we will consider the case where 
n > 2. 

It follows from (5.2) that 



Zip, x) dx = - for all x, 
p 



(5.15) 



so that 



and 



dx = 



dw(p,x) f d 



dxj 



dx 



-Z{p,x-i) [[w (0-w (x )] di 



(5.16) 
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where x° G R n is an arbitrary fixed point. 

Let us divide the domain of integration in (5.16) into two subdomains 

Q l = G r» : \£-x°\ > 5}, Q 2 = R n \Q 1 , 

where 5 > 0. This implies the decomposition (with the obvious notation) of the left-hand side 
of (5.16) into the sum v i(p, x) + v 2 {p, x). The function v i(p, x) may be differentiated under the 
integral, if x is in a small neighbourhood of x° (then \x — £| remains separated from zero). After 
that, we set x = x°, so that 



dvx(p, x°) f d 



dxi 



dx 2 



Z(p,x°-0 [[w (0 - w (x )} d£. 



(5.17) 



Let < d < 5/2, d = (0, . . . , d, . . . , 0) (d is in the j-th place). Then 

d 2 Z{p,x°-0 



v 2 (p,x° + d) — v 2 (p,x° 

dZ(p,x° + d-0 



n 2 



dx 2 



■[w (Z)-w {x )]dZ 



-\ i 



dxi 



[ Wo (Z)- Wo ( x °)]dt-± J 



dxj 



[w (0-w (x°)]d^ 



\x°-£\<2d 



\x°-£\<2d 



d 2 Z(p,x°-£) 



dx 2 



[w (Z)-w (x )]dZ 



+ 



1 

d 



2d<|a;0-^|<5 



\x°-£\<2d 

dZ(p,x° + d-Z) dZ(p,x°-0 



dxi 



dxj 



d 2 Z{p, x°-0 



dx 2 



[w (0-w (x°)]d^ 



dcf 



X 1 + X 2 + X 3 + x 4 . 



We have 



l*i| < 



C 



x° + d-£ 



-n+l+7 C 



y + d 



-n+l+7 



dy, 



\x°-£\<2d 



■\<2d 



and after the change y = d ■ rj we find that \Xi\ < Cd? — > 0, as d — > 0. Similarly, X 2 — > 0, as 
d ->■ 0. For X 3 , we get from (5.12) that 



IX, I < C 



J |x°-er n+7 ^=c 1( F^o. 



\x°-£\<2d 

In X 4 , by the Taylor formula, the expression in braces equals 

d d 3 



2dx 3 - 



Z(p,x'-£), x' = x° + 6d, 0<6< 1, 
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where 



so that 



\x - f | > |f - x°\ - \x' - x°\ > If - x°\ - d > ~|f - x°\ 



X A <Cd J \x° - f | n+7 1 d£ = Cd" 1 J \z\- n ~ 1+1 dz -> 0, 

2<2<|x°-§|<<5 2<|zj<<5(i- 1 



as d — )■ 0. 
Hence, 



dv2(p,x°) f d 2 



dxj J dx 2 

n 2 



Z(p,x°-0 [K(0 - w (x )} dC 



Taking into account (5.17) we find that, for any p > 0, x G M™, 

A x w(p,x) = J A, [Z(p,x- 0] K(0 - ^o(^)] df. (5.18) 

For each p > 0, the function Z(p,x) is a fundamental solution of the equation — Au + 
p}C(p)u = (it coincides, up to an easy change of variables, with the well-known fundamental 
solution for the equation —Am + u = 0; see Chapter 8 in |24j). Therefore AZ(p,x — f) = 
p)C(p)Z(p, x — f ), x 7^ f . Substituting this in (5.18) we find that 

Aw(p, x) = pK.(p)w{p, x) — pJC(p)wo(x) / Z(p,£)d£, 



and it follows from (5.15) that w(p, x) satisfies (5.1). ■ 

Note that a probabilistic representation of a fundamental solution of the problem (1.5) (and 
more general problems, with Levy generators instead of the Laplacian) was found in [22, 23j. 
The fundamental solutions were understood there as those of the equations obtained by applying 
the Laplace transform in time and the Fourier transform in spatial variables. For our situation, 
we obtained solutions, strong with respect to the variable x. To obtain classical solutions, one 
needs stronger assumptions regarding the function K{jp). For example, it would be sufficient to 
assume its asymptotic properties found in [H] for the distributed order case (1.3). Then the 
investigation of strong solutions would repeat the reasoning from |14j . 

Our uniqueness result for the problem (1.5) holds under much more general assumptions 
and is an immediate consequence of a deep result by E. E. Shnol (see Theorem 2.9 in [5]). Note 
that the notion of a LT-solution makes sense also for polynomially bounded solutions, that is 
such solutions w(t, x) that \w(t,x)\ < P(\x\) where P is some polynomial independent of t. 
Instead of (*), we make the following weaker assumption: 

(**) The function k is nonnegative, locally integrable, nonzero on a set of positive measure, 
and its Laplace transform JC{p) exists for all p > 0. 
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Theorem 4. Let (**) hold, and suppose that w(t,x) is a polynomially bounded LT-solution of 
the problem (1.5) with Wq{x) = 0. Then w(t,x) =0. 

Proof. The Laplace transform w(p, x) satisfies, for each p > 0, the equation Aw(p, x) = 
pJC(p)w(p,x). Thus w(p,x) is a generalized eigenfunction of the operator —A on L 2 (IR n ) with 
the eigenvalue —pJC(p) < 0. By Shnol's theorem, a nonzero polynomially bounded generalized 
eigenfunction is possible only if the eigenvalue belongs to the spectrum of — A equal to [0, oo). 
Therefore w(p, x) = 0, so that w(t, x) = 0. ■ 

Theorem 4 can be extended to some equations with coefficients depending on x, for which 
Shnol's theorem can be applied; see [51 \T2\. 
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